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Abstract 

A simple method for obtaining MCSCF orbitals and Cl natural orbitals 
adapted to degenerate point groups, with full symmetry and equivalence restric- 
tions, is described. Among several advantages accruing from this method are the 
ability to perform atomic SCF calculations on states for which the SCF energy ex- 
pression cannot be written in terms of Coulomb and exchange integrals over real 
orbitals, and the generation of symmetry-adapted atomic natural orbitals for use in 
a recently proposed method for basis set contraction. 
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I. Introduction 

For an atom or molecule with point group £, a set of functions (such as atomic 
or molecular orbitals) is said to be symmetry adapted if each function is a basis func- 
tion for an irreducible representation of Q, and the set is said to display equivalence 
if a full set of partner functions is available for each basis function. In general, for an 
atomic or molecular state that transforms according to a degenerate representation 
of £, an unrestricted SCF or MCSCF calculation on one component of the state 
will not yield orbitals which are symmetry adapted or equivalent. For example, 
SCF optimization of the fluorine atom configuration 

ls 2 2s 2 2p 2 2p 7 2p z 

will produce a 2 p z orbital with a different radial function (inequivalent) from 2 p x 
and 2p y . In addition, if a spin unrestricted approach is used, the orbitals associated 
with a or /? spin will also differ. An MCSCF optimization of a 7 D state using the 
two configurations s 2 d and sd 7 will produce optimum orbitals which are of neither 
pure s nor pure d type. Constraining the calculation so that orbitals with all the 
desired symmetry properties are obtained is referred to as imposing symmetry and 
equivalence restrictions. 

Roothaan and Bagus [l] implicitly imposed symmetry and equivalence restric- 
tions [2] in their atomic SCF calculations by solving for only the radial part of the 
orbital; the angular part was determined by 0(3) symmetry. Even in calculations 
where each orbital is expanded in the full basis set, the simple structure of SCF 
energy expressions (or, equivalently, of SCF reduced density matrices) allows the 
restrictions to be imposed via the vector coupling coefficients used to average degen- 
erate configurations- see Jackels and Davidson [3]. For example, in the F atom 2 P 
case discussed above, an energy expression corresponding to an average of energy 
associated with the configurations having the 2 p x , 2 p y and 2 p z orbitals successively 
singly occupied would be used. Jackels and Davidson also noted that if the a and 
0 orbitals were averaged, symmetry and equivalence could be imposed in a spin un- 
restricted approach. In this way many atomic and molecular states can be treated 
at the SCF level using a molecular SCF program with little or no facilities for han- 
dling symmetry: any atomic term arising from occupations of the form s m p n can be 
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treated this way, for example. There are, however, states which cannot be treated 
in this way, such those arising from certain d k occupations. The terms d 2 ( 3 F), 
d 3 ( 4 F ), d 7 ( 4 F), d 8 ( 3 F), for instance, each yield average energy expressions which 
can be expressed in terms of Coulomb ([*’t and exchange ([0*10"]) integrals alone 
only when complex orbitals are used. In terms of real orbitals, the average energy 
expressions contain integrals of the general form [tj|fc/], and such expressions can- 
not be handled by conventional SCF codes. Energy expressions which contain only 
Coulomb and exchange integrals over real orbitals can be obtained by averaging 
the F states above with the P states of the same spin from the same d occupation: 
it is actually vector coupling coefficients for these P/F averages that are given by 
Poirier et al in a recent compendium [4], although it is claimed that the coupling 
coefficients are for the F states. 

More complicated energy expressions can be handled by an MCSCF program. 
However, very few MCSCF (or Cl) programs exploit more than Dih symmetry, and 
if only one component of, say, the d 2 3 F state is optimized within D 2/1 symmetry, 
the final MCSCF orbitals will be symmetry adapted within Z? 2 /t but will not display 
atomic symmetry and equivalence. (It is relatively easy to eliminate the mixing of 
the s and d orbitals, although this will not make the five d orbitals equivalent). 
Similarly, the natural orbitals from a Cl calculation in Dih on such a state would 
again fail to show symmetry and equivalence. Such problems are not, of course, 
confined to atomic calculations, but arise in any treatment of a degenerate state 
within Dih. or its subgroups. In molecular systems symmetry and equivalence re- 
strictions could be imposed in the same manner as Roothaan used in his work on 
atoms, that is, a special code for each symmetry of interest could be developed. 
While this might be consider an elegant approach, it is impractical since it would 
lead to the problem of developing and supporting a large number of very complex 
codes. In the present work, we first show formally how symmetry and equivalence 
restrictions are imposed, and then describe how these can be implemented compu- 
tationally. A simple example is given to illustrate the method. This approach to 
symmetry and equivalence restrictions is a very powerful one, and its implementa- 
tion is a much simpler task than that of extending codes designed originally for Dzh 
and its subgroups to include higher symmetries. The ability to impose symmetry 
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and equivalence restrictions is of particular current importance in the construction 
of new atomic natural orbital (ANO) basis sets. 

II. Symmetry Adapted Natural Orbitals and Geminals 

The symmetry properties of density matrices and natural orbitals have been 
reviewed extensively by McWeeny and Kutzelnigg [5] and by Davidson [6], and we 
shall generally follow their treatment. The kernel of the electron density opera- 
tor (t denotes both spin and spatial coordinates) is given in terms of an 

orthonormal one-electron basis {0} as 

( 1 ) 

pq 

where ~i pq is an element of the density matrix. This is a Hermitian matrix and 
can be brought to diagonal form by a unitary transformation, yielding natural spin 
orbitals (NSOs) {x}: 

'Kg r') = n p^( r )Xp(’’ / )» (2) 

p 

where n p is the occupation number of natural spin orbital Xp- The NSOs are of 
pure ct or j3 spin type if the Hamiltonian contains no spin-dependent operators [5]. 
That is, a given NSO can be written as a product of a spatial orbital 4> p (r) (r 
denotes spatial coordinates) and a spin factor a or 0. However, it is not generally 
the case that a given 4> p (r) will appear with both a and 0 spin factors: the partner 
function (within the given irreducible representation of the spin group) of an NSO 
is generally not itself an NSO. This is, in effect, a loss of equivalence properties 
under the spin group for the NSO. Only in the case of zero total spin projection 
(Ms = 0) is equivalence obtained. This is obviously a nuisance, as different NSOs 
will be obtained from density matrices for different Ms values. It is therefore 
customary [5,6] to define a spin-free density matrix (denoted p) as the sum of the 
a-spin and /?-spin blocks of 7. It is this matrix which is normally computed as 
the “density matrix” in electronic structure codes, and its eigenvectors are termed 
natural orbitals (NOs). As will be seen, this spin- averaging is exactly analogous to 
imposing spatial symmetry and equivalence restrictions. 
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The kernel of the spin-free density operator can be written as 

/>(>•;»•') = ( r T- (3) 

It is convenient to assume that the set of orthonormal orbitals {<p} possess the de- 
sired symmetry and equivalence properties. (This is not a practical restriction, as 
some such set of orthonormal orbitals can always be generated straightforwardly, 
say, by diagonalizing the one-electron part of the Hamiltonian (h). See for example 
Ref. 7, where it is shown that it is also trivial to sort the eigenvectors of h by irre- 
ducible representations). We denote by <p$ a a basis function for row a of irreducible 
representation p\ p is simply a counting index. Equation (3) thus becomes 

p( r > r ') = £EX>rM<rxvr. <»> 

fiu ab pq 

where p ao, ^ b is a sub-block of p. What is now desired is to obtain NOs with full 
symmetry and equivalence properties from p. Unfortunately, the eigenvectors of p 
will only display these properties when the wave function transforms according to a 
non-degenerate irreducible representation of £ [5,6]. In other cases, it is necessary 
to project out p°, that component of p that transforms according to the totally 
symmetric irreducible representation of Q, and use the eigenvectors of p as the 
NOs. These orbitals display full symmetry and equivalence properties and differ 
as little as possible (in a least-squares sense [5]) from the eigenvectors of p itself. 
The following projection approach, which is essentially that used by Davidson [6], 
is easily implemented computationally. 

The component of p which transforms as a basis function for row a of irreducible 
representation p is obtained [6,8] as 

P l “‘=g-'n lt 'E D Zt.( R Y R l’ Ir '’ Re S’ ( 4 ) 

R 

where g is the order of § and D M (i?) is one of a set of unitary irreducible repre- 
sentation matrices for p, which is of dimension n M . A set of basis functions for all 
rows of irreducible representation p is obtained by choosing a linearly independent 
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subset of the functions obtained using all possible values of index 6 in (4). For the 
totally symmetric irreducible representation we have 

f° = g-'Y. R P R ~'- (5) 

R 

It is simplest to project p block by block: for block (p 0 ) Ma,I/b we obtain 

(/)<•«.-* = = g -'YY.<> l ‘ c '‘"‘ Dl iA R )DU R Y ■ ( 6 ) 

R R cd 

But 

9 ~‘ £DSA R )WA R y = n; ( 7 ) 

R 

by the great orthogonality theorem [8]. Substituting ( 7 ) in (6) gives 

1 Y, p"' 1 "' (8 ) 

C 

for the only non- vanishing blocks of p°. This is very simple to implement computa- 
tionally: after identification of blocks of the density matrix by their transformation 
properties under the diagonal blocks are averaged within representations, while 
all off-diagonal blocks are discarded. p° can then be diagonalized, to give NO oc- 
cupation numbers n p and NOs <j>£ a : 

?°(r;r') = EEEO'XI'T. (9) 

P a p 

It is easily seen that projection of that component of -7 that transforms accord- 
ing to the totally symmetric irreducible representation of the spin group (together 
with a renormalization so that NO occupation numbers He between zero and two) 
corresponds exactly to the spin averaging procedure described above. 

The above approach is sufficient to allow symmetry and equivalence restricted 
NOs to be generated from a given density matrix. However, in a general energy ex- 
pression the second-order reduced density matrix will also appear, and it is therefore 
desirable, for the purposes of averaging energy expressions, to extend the above ap- 
proach to the second-order reduced density matrix. The kernel of the second-order 
reduced density operator [ 5 ] is 

r(ri,T 2 ;rj',r2) = V>m(»i )^ p {t 2 )T mn Pq ^n{r[) ‘^*(*2) *• (10) 

mnpq 
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Here, as in (l) above, {t/>} is a set of orthonormal spin-orbitals, and T is the second- 
order reduced density matrix. Note that the indices on T have been ordered to cor- 
respond to those on two-electron integrals given in charge density notation. Again, 
in order to avoid obtaining different results for different Ms values it is usual to ob- 
tain a spin-free kernel P(ri,r 2 ;r'j,r 2 ) by averaging over spin blocks. Full details are 
given by McWeeny and Kutzelnigg: P and the associated reduced density matrix 
P 

P(ri,r 2 ;ri,r^) = ^m(ri)v?p(r 2 )P m np g V?n( r i) V«(»a)*, (H) 

mnpq 

are the quantities computed in density matrix-driven MCSCF codes. 

In terms of the orthonormal symmetry-adapted basis {<p£°} introduced above 
equation (11) becomes 

/ > (r„r 2 ;r'„r;)= £ £ £ *>K‘(flX(r J )/>^* cW *>S*(r' 1 )V«(rJ)*. (12) 

abed mnpq 

P°, the component of P that transforms according to the totally symmetric irre- 
ducible representation, is obtained by projection of symmetry blocks of P. For a 
given block we have 

g~l ^ RRPP aUbKcXd R~ } R~* = 0 _1 ^ plLivjKkXl 
R ijkl 

* (Y KiWWkW D^RY DlARY). 

R 

The factor in braces in (13) can be regarded as a product of elements of unitary 
reducible representations for the direct products /x 0 k, and u ® A. The unitarity of 
these reducible representations gives 

Y D^(R)D^(R)Dii(RY Di , (*)• = ■*«».«. ( 14 ) 

R 

Insertion of (14) into (13) yields the final symmetrization formula, however, it is 
obvious from inspection of (14) that the final form will be more complicated than for 
the first-order case (8). In the following section we will show how the symmetrization 
of P can be avoided in MCSCF calculations. If natural geminals with symmetry 
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and equivalence restrictions were required [5], this step would have to be performed, 
of course. 

III. Illustration 

As a simple example of using equation ( 8 ) to obtain symmetry and equivalence 
restricted NOs from a Cl wave function, we consider the excited 2 E' state of planar 
CH 3 . The molecule has Dzh symmetry with SCF occupation 

la' 2 2a' 2 la' 2 ' 2 le' 3 . (15) 


In a program which handles only D 2 /i and its subgroups, this system must be treated 
using C? v symmetry. The two components of the 2 E' state (15) in Cz v symmetry 
are 

2 A\ Ia 2 2a 2 l6 2 3a]l6 2 (16) 


and 


2 Bi Ia 2 2a 2 l6 2 3a 2 l6}. 


(17) 


First, we note that if the orbitals have full D^h symmetry, a Cl calculation which 
includes all single and double excitations from both (16) and (17) will yield the same 
energy for the 2 A\ and 2 B\ states. The first-order reduced density matrices from 
these two Cl wave functions are given in Table 1 . The two problems which prevent 
the eigenvectors of either of these density matrices displaying Dzh symmetry are 
immediately apparent. First, density matrix elements such as pmn 1 are not equal 
to pmn’ in either matrix. Second, both matrices contain non-zero elements of the 
form pmn 1 • It is straightforward to see that the projection procedure defined by 
equation ( 8 ) above eliminates both problems, and the same projected density matrix 
is obtained from both of the original matrices. At this stage it should be pointed 
out that in practical calculations it is important to ensure that phase relationships 
among degenerate orbitals are consistent. That is, applying the shift operator P£ b , 
which generates a basis function for row a of irreducible representation p from a 
basis function for row 6 , must yield 



( 18 ) 
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where A is the same phase factor, for all values of the counting index p. Although 
this condition does not affect the computed energy, it is important in averaging to 
obtain symmetry and equivalence properties. 

It can also be seen from Table I that the same averaged density matrix as ob- 
tained by projection can be obtained by averaging the density matrices for the two 
components of the 2 E' state, since the two wave functions and the corresponding 
density matrices are related by a shift operator. However, if the averaged density 
matrix is obtained from averaging calculations on the individual components, it is 
of paramount importance in practice that the component wave functions be related 
via a shift operator. For example, it was specified above that the CH 3 wave func- 
tions comprised all single and double excitations out of both components, (16) and 
(17). This condition is sufficient to guarantee that the Cl wave functions will be 
related via a shift operator, even though it means that single and double excitations 
out of reference configurations of the “wrong” symmetry must be included. If such 
excitations are excluded, then the wave functions for the two components will not 
be related by a shift operator. This makes it impossible, for example, to impose 
first-order interacting space restrictions [9] on the Cl wave functions. However, 
the symmetrization procedure of equation (8) requires only one Cl calculation to 
be performed, and first-order interacting space restrictions can be imposed, as the 
symmetrization procedure effectively generates the result of applying a shift opera- 
tor to this Cl wave function. It therefore seems preferable to symmetrize a density 
matrix for a single component of a degenerate state rather than to average density 
matrices for all components. 

The ability to obtain Cl NOs with full symmetry and equivalence properties is 
very important for generating contracted basis sets for molecular calculations using 
atomic natural orbitals (ANOs) and a general contraction scheme [10]. Because 
of the shell structure exploited in efficient integral codes [ll], it is vital that the 
contracted basis display full atomic symmetry and equivalence, and it is not always 
possible to find a suitable non-degenerate state of a given atom for generating such 
ANOs directly. 

Our computational implementation of equation (8) for obtaining symmetry- 
adapted NOs is very simple and we briefly outline it here. We refer to the full 
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symmetry it is desired to impose as “high symmetry” and the appropriate subgroup 
of Dih as “low symmetry”. Required as input is the number of high symmetry 
irreducible representations, their degeneracy, and the number of orbitals in each 
irreducible representation. Then, for each row of each irreducible representation 
(na ) , the sequence numbers of the low symmetry orbitals which transform as basis 
functions for /xa are given. The program loops over distinct pairs of irreducible 
representations p and v. If p and u are not equal, />£“•*'* is set to zero. In the case 
of ii = is, ppg ,flb is set to zero for all a ^ b, while for a = b, we average 
over all o. Diagonalization of p produces NOs with full symmetry and equivalence 
restrictions. This implementation only requires the identification of the equivalent 
(high symmetry) orbitals in the original orbital list. 

Imposing symmetry and equivalence and restrictions in the NOs in this manner 
is similar to the method by used in the SCF Fock matrices, where the equivalent 
blocks are averaged and the off-diagonal elements set to zero. Also this method is 
being used in MCSCF treatments based on NOs and the BLB theorem [12]. 

IV. Symmetry and Equivalence Restrictions in MCSCF Calculations 

We now turn to the question of imposing symmetry and equivalence restrictions 
in second-order MCSCF calculations. In most MCSCF methods [13] the variation 
of the MCSCF energy with respect to orbital rotations is written in terms of first 
and second-order reduced density matrices, and the methods of section II could 
therefore be used to impose symmetry and equivalence restrictions on the MC- 
SCF optimization. It is possible to avoid the more complicated symmetrization 
of the second-order reduced density matrix, however. The MCSCF energy can be 
expanded in the orbital rotations as 

E = Eo + QpqXpg -I- - ^ Hp qt r e Xp q X re , (19) 

P>q p>q r>e 

where Eq is the energy with the current MOs, g is the gradient vector, H is the 
Hessian matrix and the antisymmetric matrix X is used to parametrize the orbital 
rotations. Seeking a stationary point of (19) leads to the Newton-Raphson equations 

^ Hpq, T sXrt ~ ~9pq- ( 20 ) 

r>t 
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Symmetry and equivalence properties will be maintained only if the rotations are 
restricted to mix orbitals of the same symmetry type, (that is, if the only non- 
vanishing X rs are those for which r and 5 transform as basis functions for the same 
row of the same irreducible representation), so this is the condition that must be 
imposed on the Newton-Raphson equations (20). Labelling the symmetry blocks in 
(20) explicitly, we have 


\ ^ V ' \ ^ 11 tiana,vbi/b \ri/bvb 

~ ~ 9 P<1 • 


(21) 


V b T>8 

The symmetry structure of H is clearly much simpler than that of the second-order 
reduced density matrix in equation (12), as the blocks which transform according 
to four different symmetry indices do not appear in (21). 

For the MCSCF case, our implementation is similar in spirit to that for the 
first-order density matrix; we wish to specify only the number of high symmetry 
irreducible representations and, for each irreducible representation, its degeneracy 
and the list of the low symmetry orbitals associated with each of its rows f. 1a . It is 
obvious that the gradient vector can be symmetrized using exactly the same scheme 
as for the first-order reduced density matrix. For the Hessian, it is desirable first 
to compress the notation somewhat. As we are interested only in rotations within 
symmetries, we can use a single label fj.a, and we can represent distinct pairs rs 
by a compound counting index t. The blocks of X of interest are then those with 
elements Xf°. Given the high symmetry irreducible representation and orbital 
list information, potential orbital mixings which would correspond to symmetry 
breaking can be identified. It is now possible to generate a list of allowed mixings 
for each irreducible representation, and to use this rather than the list information 
about the original orbitals. The elements of the Hessian which correspond to a 
mixing which breaks symmetry are now set to zero (excluding diagonal elements 
which can be set to an arbitrary positive value). Note that the elements of g 
corresponding to rotations which break the desired symmetry are eliminated in 
the projection of the gradient. The program loops over distinct pairs of irreducible 
representations \i and v. For the case of n — */, and a = 6, the Hessian blocks 
are averaged over all a. This is equivalent to the processing of the gradient or first- 
order reduced density matrix with the orbital indices replaced by the compound 
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mixing indices. For a ^ b, the elements H? a ^ b + Ht^ a ' t ^ b are averaged for each value 
of t > w over all a > b. Thus for an atom with d orbitals, the ten lower-triangular 
off-diagonal d-d blocks and their transposes would be averaged. For /z ^ v, the 
elements H^ l/b are averaged for all a and b. Thus for an atom, fifteen equivalent 
p — d blocks would have to be averaged. 

In the our implementation, we explicitly store the equivalent blocks of the Hes- 
sian, thus retaining its orginal low symmetry dimension. Clearly the dimension of 
the Hessian could be reduced to the unique orbital mixings in the high symmetry. 
However, by retaining the full low symmetry Hessian dimension the program modi- 
fications are limited to the symmetrization of the Hessian. The described method to 
maintain orbital symmetry and equivalence during MCSCF optimizations is similar 
to that used by Ruedenberg, Cheung and Elbert in the ALIS program [12]. 

With the same techniques it is also straightforward to impose symmetry and 
equivalence restrictions in MCSCF calculations which include orbital-CI coupling: 
the Cl gradient and Cl Hessian can be constructed using one component of a de- 
generate state, while the Cl/orbital coupling Hessian can be symmetrized in the 
same manner as the orbital gradient. 

This approach can be used to obtain SCF solutions for the d occupations dis- 
cussed in the introduction. For example, for d 3 there are ten determinants with 
S = Ms — 3/2, one 4 A g and three 4 B i 9 , 4 i?2 S , and 4 £ 3 g in Dzh symmetry. Op- 
timum SCF orbitals for the 4 F atomic state can be obtained by performing an 
MCSCF calculation with atomic symmetry and equivalence restrictions imposed, 
solving for any of the seven possible components: 4 A g or the two lowest (degenerate) 
4 B i g , 4 B 2 g, or 4 Bz g solutions. All seven solutions will be exactly degenerate and 
will yield identical orbitals. The third root in each of the B symmetries corresponds 
to a component of the d 3 ( 4 P) state. 

It is well known [14] that the orbitals for the lowest state arising from the 
3d n 4s 2 and 3d n+1 4s 1 occupations in transition metals are very different. Using 
the approach described above, we have recently performed state-averaged CASSCF 
calculations with symmetry and equivalence restrictions for the 3 F (3d 2 4s 2 ) and 5 F 
(3d 3 4s J ) states of Ti. These CASSCF calculations were followed by Cl calculations, 
and the Cl first-order density matrices for the two states were averaged. Since full 
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symmetry and equivalence restrictions were imposed, these averaged NOs can be 
used to determine an ANO contraction. As discussed by Shavitt [15] in his review 
of configuration interaction methods, averaged NOs are expected to be a good way 
of obtaining an equivalent treatment of two states, and for Ti it is found that the 
average NOs yield a more equivalent description of 3 F — 5 F separation than either 
3 F or S F orbitals. A [5s 4p 3 d If] contracted set yields a 3 F — 5 F separation 
which agrees to within 0.03 eV with the result obtained using the uncontracted 
(14s lip 6d 4f) basis set. One example of the need for methods for performing 
beyond-SCF calculations in which symmetry and equivalence restrictions are im- 
posed is in the construction of such compact ANO contractions. 

While the ANO contractions are one clear example of the need to impose 
symmetry and equivalence restrictions, it is not always an advantage to use the full 
molecular symmetry. One example would be the distortion away from an isolated 
point of high symmetry where the wave function is degenerate. The imposition 
of symmetry and equivalence at the high symmetry point, but not at the lower 
symmetry geometries could result in a discontinueous potential surface — see the 
discussion by Davidson and Borden [16]. One solution is to state-average for the two 
states, then at the high symmetry point, as discussed above, the density matrix has 
the full high symmetry, and the potential energy surfaces are smooth. While we have 
sucessfully used this approach in one application, the question of symmetry breaking 
is very complex and it is not yet clear how general a solution state-averaging will 
be. 

V. Conclusions 

We have presented a simple method of modifying a Z> 2 fc-based MCSCF/CI 
program to obtain full symmetry and equivalence restricted MCSCF solutions and 
Cl natural orbitals. Among other advantages, this allows ANO basis sets to be 
generated using standard program systems, and SCF solutions to be obtained for 
cases where conventional SCF approaches cannot be used. 
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Table I. The MO density matrix for the 2 E' X and 2 E' y components of the 2 E' state 
state of CH3. The density matrix is blocked and labelled by Ci v , the symmetry 
used in the Cl calculation, while the orbitals are labelled using £> 3 ^, the symmetry 
of the molecule, and that used in the SCF calculation. 


«i a 

1.987957 

1.987957 


Ol 



e'x 

-0.003925 

1.002602 





0.003925 

1.988767 




a'i 

0.003455 

0.030386 

0.010319 




0.003456 

-0.030386 

0.010319 



a\ 

0.003404 

-0.018831 

-0.000721 

0.005609 



0.003404 

0.018831 

-0.000721 

0.005609 


c x 

-0.009404 

0.031865 

0.000381 

-0.001446 

0.006198 


0.009404 

-0.043940 

-0.000381 

0.001446 

0.005489 

e v 

1.988767 

1.002602 


b 1 



e 'v 

-0.043940 

0.005489 





0.031865 

0.006198 




a'i 

1.972070 

1.972070 


b 2 



a” 

-0.011975 

0.020990 





-0.011975 

0.020990 



- 


a The top line is for the 2 E' X component while the bottom line is for the 2 E' y 
component. 
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